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Работа посвящена описанию и экспериментальному исследованию возможностей алгоритма 
обработки слабоконтрастных изображений с целью улучшения их яркостных характеристик, а 
именно, повышения контраста для обеспечения достоверности визуального анализа. Цель достигается 
за счет использования аппарата нечетких множеств, применения локальных и глобальных 
статистических характеристик при использовании оконных преобразований. Представлены 
результаты экспериментальных исследований возможностей предложенного алгоритма на примере 
реальных медицинских изображений. 
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Введение 

Изображения, сформированные различными информационными, следящими, 
диагностическими системами очень часто имеют качество, недостаточное для 
проведения достоверного анализа. Недостаток качества изображения обусловлен 
различными причинами, такими как шум, дефекты систем формирования 
изображений или особенностями объекта исследования. 

Улучшение качества изображения — важный, а иногда и судьбоносный, как в 
случае медицинских исследований, фактор в принятии решения. Сложность этой 
процедуры обусловлена как существенным различием характеристик анализируемых 
изображений, так и разнообразием поставленных целей. Кроме того, отсутствие 
объективных критериев оценок качества исходных данных и конечного результата 
приводит к сложности выбора алгоритмов для решения этой задачи [1]. 

Постановка задачи 

Метод усиления контрастности с использованием функции протяженности 
гистограммы эффективно используется в обработке широкого класса изображений. 
Учитывая характеристики окрестностей, удается идентифицировать участки 
изображения по уровню контрастности и соответствующим образом их 
обрабатывать. Благодаря этому достигается более тонкая обработка мелких деталей. 
Однако изображения должны отвечать двум требованиям. Они не должны содержать 
большого количества импульсных выбросов и темные или светлые участки большой 
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площади. В первом случае это может привести к неадекватному вычислению 
функции протяженности гистограммы, а во втором — к неэффективному усилению 
контрастности. Поэтому, если изображение не отвечает указанным выше 
требованиям, следует провести его фильтрацию или (и) градационную коррекцию. 

Классическая обработка контрастности обычно основывается на глобально 
определенной функции протяженности. Традиционно выравнивание гистограммы — 
это тоже глобальная методика в том смысле, что усиление основывается на 
выравнивании гистограммы всего изображения. Однако, хорошо известно, что 
использование только глобальной информации часто недостаточно для достижения 
хорошего усиления контрастности (например, глобальные методы вызывают эффект 
интенсивной насыщенности). 

При  адаптивном усилении контрастности новая — интенсивность 
устанавливается для каждого пикселя в соответствии с адаптивной передаточной 
функцией, разработанной на основе локальной статистики (значений локального 
минимума/максимума/среднего по интенсивности) [2]. 

Анализ последних достижений и публикаций 

Нечеткая обработка изображений является своего рода нелинейным 
преобразованием, отличие которого от других хорошо известных методик 
заключается в том, что оно выполняется над функциями принадлежности пикселей к 
предварительно определенным кластерам [3]. 

Таким образом, первым этапом обработки является фаззификация — генерация 
подходящих значений функции принадлежности. 

Этап фаззификации может интерпретироваться как тип кодирования входных 
данных, который необходим еще и потому, что исходные значения часто являются 
неполными и/или неточными. Кроме случайности, которая может управляться в 
соответствии с теорией вероятности, есть три других вида неопределенности в 
задачах обработки изображения: двусмысленность серости, геометрическая 
нечеткость, неопределенность знания. 

Как и в других областях применения теории нечетких множеств [4] шаг 
фаззификации во многом определяет конечный результат. Число, форма и 
расположение каждой функции принадлежности может/должна быть адаптирована к 
конкретной задаче. Возможны следующие варианты фаззификации: 

— локальные — результирующее значение зависит лишь от значения в 
соответствующем входном пикселе изображения; 

— фокальные (оконные) — результат определяется малой окрестностью каждого 
пикселя (они играют особую роль в силу пространственной автокорреляции 
пикселей); 

— зональные (глобальные) — при преобразованиях используются интегральные 
характеристики изображения. 

Часто используемая фаззификация на основе анализа гистограммы относится к 
глобальным методам. Данный подход предполагают учет уровня серого каждого 
пикселя, на основе чего происходит определение функции принадлежности к одному 
или несколько классам, например, таких как очень темный, слегка яркий, средний и 
тд. Для его использования необходимы некоторые априорные знания об 
анализируемом изображении, такие как минимум и максимум частот уровней 
серого, среднее значение. 

Методы сегментации, фильтрации шума и т.п. работают на предопределенной 
окрестности пикселей. Чтобы использовать нечеткие подходы к таким операциям, 
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этап фаззификации также должен быть выполнен в пределах выбранной 
окрестности. Такой подход требует больше времени вычислений по сравнению с 
применением гистограмм. Кроме этого, в таком случае большее влияние оказывают 
шумы и выбросы, которые могуг привести к ошибочным значениям функции 
принадлежности. 

Как показали исследования [5], пространственная корреляции пикселей, как 
правило, в достаточной степени учитывается в окрестности размером 3х3. Учет 
соседства может осуществляться различным образом в зависимости от поставленной 
задачи, как с учетом комбинаций локальных контрастов отдельных элементов, когда 
все пиксели считаются равнозначными, или окрестностей (ближайших соседей) 
пикселей, так и на базе интегральных характеристик всего изображения — полноты 
использования уровней серого, гистограммы, скорости изменения яркостей и т.д. 

Сформированные значения функций принадлежности модифицируются с 
помощью соответствующего нечеткого подхода (например, на основе использования 
усиливающего оператора [6]), результатом чего является новые значения функций 
принадлежности каждого пикселя, на основе которых синтезируется 
результирующее изображение [7]. 

Цель исследования 

В данной работе предложен метод адаптивного преобразования 
слабоконтрастных изображений, который позволяет обеспечить улучшение их 
яркостных характеристик, а именно, усиления контраста, с целью повышения 
достоверности визуального анализа изображений при принятии решений 
диагностики, за счет использования аппарата нечетких множеств, локальных и 
глобальных статистических характеристик в процессе выполнения оконной обработки. 

Изложение основного материала 

Представленный в этой статье алгоритм адаптивной обработки изображений 
представляет собой модификацию (на основе применения фокальных и глобальных 
статистических показателей) метода, основанного на применении усиливающего 
оператора [6], и состоит из следующих шагов: 

1. Масштабирующее преобразование каждого цветового канала исходного 
изображения на [0,1]. 


2. Каждый цветовой канал исходного изображения разбивается на 
неперекрывающиеся окна заданного размера. Выбор размера окна зависит от 
исходных данных и целей анализа. Для большинства изображений рекомендуется 
использовать размеры окон от 3х3 пикселей до 6х6 пикселей. 


3. Для пикселей текущего окна формируется значение 4х, , по следующей формуле: 


Е 
пах 
У = 
‚У 
Иху = 1 : (1) 
х, у Е) 
где их, - пиксели текущего окна, а и ‚ и Е. вычисляются по 
следующим формулам: 
итах — [утих +. улих /о, (2) 
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и (и т }/2 


Е 

Ер =2+и +[-и, (4) 
причем у”, и", ри Г, Г, 1 - максимальные, минимальные и 
средние значения для текущего окна и выбранного цветового канала исходного 


' 


изображения, соответственно, а вычисляется по формуле: 


и =(+7+0.5)/3. (5) 
1. Осуществляется пропорциональное повышение яркости и контраста на 
основании применения следующего степенного преобразования: 


Ра = (3) 


а = АЕ Их м . (6) 


2. Вычисляются значения функции принадлежности Й у по формуле: 


1 —и+0.5 1 
2 0 < д», у < 0.5 


х,у 
— _ 7 
1 +и-0.5 ий 


Их,у = 
1-2 дна! 

Данное преобразование, по сути, осуществляет повышение контраста за счет 
пропорционального сдвига пикселей, интенсивность которых не превышает 0.5, к 0, 
а остальных — к 1. Однако использование такого подхода также приводит к 
снижению яркости для пикселей, интенсивность которых не превышает 0.5, что и 
обуславливает необходимость выполнения предыдущего шага. 


3. Происходит формирование пикселей текущего окна выбранного цветового 


1 
канала выходного изображения ух, на основании следующего преобразования: 


—— —1 
и =иах — К, РУС -1 |. (8) 


Представленный модифицированный алгоритм был применен при обработке 
различных полутоновых медицинских изображений, примером которых является 
снимок, представленный на рис. 1. 

На рис. | а приведено полутоновое медицинское изображение — томограмма 
головного мозга, выполненная с целью диагностики наличия гематомы, а также 
определения области ее влияния в случае выявления. Гистограмма исходного 
изображения (рис. 1 6) показывает, что изображение нельзя отнести к классу 
низкоконтрастных. Однако диагностирование гематомы по исходному снимку 
затруднено из-за ее расположения на низкоконтрастном участке изображения 
(фрагмент, обведенный прямоугольником), а область ее влияния вообще визуально 
неразличима. Таким образом, это изображение является слабоконтрастным. 
Введение рентгеноконтрастного вещества (рис. 1 в) позволяет более четко выделить гематому. 
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а | б В 


Рис. | Рентгеновская томограмма головного мозга: а — исходный полутоновой снимок 
(204х201); б — его гистограмма; в — результаты введения рентгеноконтрастного вещества. 


На рис. 2 представлены результаты обработки снимка (рис. 1 а) исходным и 
модифицированным методами. Применение исходного алгоритма с размером окна 4х4 
пикселей (рис. 2 а) приводит к существенному ухудшению визуального анализа из-за 
размытия объектов в области интереса. Использование же окна, равного размеру всего 
изображения (рис. 2 6) хоть и выделяет гематому, но не позволяет выделить область ее 
влияния. Применение предложенного модифицированного алгоритма с окном 4х4 
позволяет как выделить гематому, так и область ее влияния (рис. 2 в). Хотя следует 
заметить, что недостаточный уровень яркости в области интереса не позволяет четко 
определить границы области влияния гематомы. 

Использование окон небольшого размера в предложенном методе обусловлено 
также борьбой с «оконным эффектом» — различимости стыков соседних окон, что может 
затруднять процесс визуального анализа полученных результатов. Использование же 
слишком больших окон (например, совпадающих с размером снимка) приводит к 
получению результирующего изображения недостаточной яркости (рис. З а). 

На рис. 3 б также показана важность применения 4 шага предложенного алгоритма 
(пропорциональное повышение яркости). Если при обработке исходного снимка 
предложенным методом с окном 4х4 данный шаг опустить, то это приводит к 
существенному снижению яркости получаемых результатов, что затрудняет визуальную 
идентификацию объектов интереса. 

Предложенный алгоритм также может быть применен при обработке цветных 
изображений и, как показали экспериментальные результаты, в случае применения окон 
небольших размеров обеспечивает сохранение цветового баланса и повышение контраста 
(ввиду черно-белой печати полученные в результате изображения не приводятся). 


а б В 


Рис. 2 Обработка исходного снимка: оригинальным методом: а — окно 4х4; 
б — окно совпадает с размером снимка; в — модифицированным методом (окно 4х4) 
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Рис. 3. Обработка исходного снимка модифицированным методом: 
а — окно совпадает с размером снимка; б — без выполнения шага 4 (окно 4х4). 


Выводы 

Предложенный в данной работе модифицированный алгоритм адаптивной 
обработки слабоконтрастных изображений позволяет добиться лучшего качества 
получаемых результатов по сравнению с исходным методом при использовании 
оконных преобразований. Предложенный метод обеспечивает повышение контраста, 
а также некоторое повышение яркости (зависит от исходного снимка) как для 
полутоновых, так и для цветных изображений. Достоинствами предложенного 
метода является адаптивность, отсутствие управляющих параметров и относительно 
невысокая вычислительная нагрузка. К недостаткам следует отнести недостаточный 
уровень яркости получаемых изображений, возможность нарушения цветового 
баланса при обработке цветных изображений, а также невозможность работы с 
темными снимками из-за особенностей процедуры повышения контраста. 
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ВЕЗОМЕ 

Г.С. Асбте& та, А.А. Уе2огоу 

АЧдарбуе те®о4 0 пиргоуе 10\-сопбга${ ппазез$ Бабед оп Фе йл72у @гапюгтпавоп 

Ть$ агасе 1$ деус{е4 ю дезсириоп ап ехрегипеша! гезеагсвез оЁ фе абШиез оЁ 
Фе адариуе тефо4 оф 1о\-сопёгаз( ппазез ргосеззте Базе оп из? Фе ЁРа2лу 5е{5. ТВе 
ригрозе о? 1$ ргосезз1те 15$ (0 епрапсе Фе Бе тез$ сВагасет$с$ оЁ фе опотпа! ипазе, 
пате[у, епВапстя сопга$ {ю паргоуе Фе гейаб у о? у1зча| апа[у$15. 

Тре Био тез$ оЁеасВ рахе! оЁ Фе сштепе ипап$Рогтайоп \1п9доу\ 1$ 5е1 Базе оп Фе 
адариуе бап$Ёег Гапсйоп. ТЬ1$ ГалпсНоп 1$ саси]меа Бу изше Фе 1оса! ап оба! 
ай$Иса| спагасет1$Яс$ Гог еасВ со]ог сваппе] (ве пииитит / фе тахипит / Фе ауегасе 
Био тез$ уащез). 

Тре ехрегипета гезеагсВез$ оЁ фе абШ@е$ оЁ Фе ргорозе4 азот детопзгае 
фе ппргоуетепе оЁ ппазе диаШбу ш сотрайзоп УИ опршта| тефо4. ТЬе ргорозеа 
а]еогифт ргоу14ез сопгаз( ппргоуетеп( ап зоте Би тез$ епрапсетепи (4ерепа$ оп 
Фе зоигсе птаге) Гог Фе со]ог ап отаузсае ппазез. Тре а4уатщасез оЁ Ше ргорозеа 
тешфшоЧ аге Фе адараБШиу, Фе абзепсе оЁ сопёо| рагатеег$ ап теануеу 1о\ 
сотршавйоп со$е. Тре гези {$ оЁ Ше ехрегитегца| гезеагсВез аге оМате4 Гог геа| пе1са1 
ппазез ргосезз1п®. 
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